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Abstract 

We discuss a class of flat FRW cosmological models based on D=4 axion-dilaton 
gravity universally coupled to cosmological background fluids. In particular, we in- 
vestigate the possibility of recurrent acceleration, which was recently shown to be 
generically realized in a wide class of axion-dilaton models, but in absence of cosmo- 
logical background fluids. We observe that, once we impose the existence of radiation 
- and matter - dominated earlier stages of cosmic evolution, the axion-dilaton dynamics 
is altered significantly with respect to the case of pure axion-dilaton gravity. During 
the matter dominated epoch the scalar fields remain either frozen, due to the large 
expansion rate, or enter a cosmological scaling regime. In both cases, oscillations of 
the effective equation of state around the acceleration boundary value are impossible. 
Models which enter an oscillatory stage in the low redshift regime, on the other hand, 
are disfavored by observations. We also comment on the viability of the axion-dilaton 
system as a candidate for dynamical dark energy. In a certain subclass of models, 
an intermediate scaling regime is succeeded by eternal acceleration. We also briefly 
discuss the issue of dependence on initial conditions. 



1 Introduction 



Although no fundamental scalar particles have been discovered yet, many attempts to 
extend the Standard Model of particle physics naturally introduce new scalar and/or 
pseudo-scalar degrees of freedom, e.g. scalar superpartners of Standard Model fields 
in a supersymmetric framework, moduli fields related to geometric properties of com- 
pactified extra dimensions, etc. Whereas on the one hand this leads to challenging 
problems [T], on the other hand these new fields could provide interesting candidates 
of dynamical dark energy [2] . 

Of particular interest in cosmology are low-energy effective theories where the gravi- 
tational sector includes, apart from the metric tensor, scalar (or pseudo-scalar) degrees 
of freedom that could provide a gravitational interpretation of early (inflationary), or 
recent (quintessential), periods of acceleration. 

In the extensively studied case of Scalar-Tensor (ST) theories of gravity [3], the 
existence of one or more scalar partners of the graviton leads to modifications of the 
Hubble expansion and of Newton's law [U [5]. While ultra-light scalar fields are in 
general potentially dangerous sources of new long range forces, interestingly enough, 
ST theories are protected against any violation of the weak equivalence principle by a 
universal metric coupling between matter and the gravity sector [5]. Because of this 
property, ST theories provide a natural framework to address the issue of dynamical 
dark energy On the other hand, inflation can also be successfully achieved in a ST 
picture [7j. 

The ST scenario can be generalized to include also pseudo-scalar fields. This is the 
case, for instance, of axion-dilaton (AD) gravity, which can be viewed as a prototype 
of theories where a "dilaton-like" scalar and an "axion-like" pseudo-scalar appear as 
spin zero partners of the graviton. Such a picture naturally emerges from the Neveu- 
Schwarz bosonic sector of the low-energy string effective action |8j. Black hole solutions 
of this theory have been found in [9] , while domain wall solutions were given in . A 
more general class of stationary supersymmetric solutions was discussed in [TT]. Cos- 
mological implications of such a theory - up to first order in perturbation theory - were 
investigated in [12], where for the first time the spectrum of cosmological perturbations 
was computed. 

It has been pointed out recently that AD gravity theories with an exponential dilaton 
potential admit cosmological solutions which give rise to the interesting phenomenon 
of recurrent acceleration [H] . By a detailed phase-space analysis of the AD dynamical 
system, the authors verified the generic occurrence of recurrent acceleration in the 
regime of a spiral focus associated to a runaway behavior of both fields. In this picture, 
the present acceleration does not appear as a peculiar stage of the cosmic history, being 

* Exhibiting a mass of the order of the present value of the Hubble parameter 
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likely a transient or even recurring phenomenon. In particular, they conclude that the 
future evolution of the universe is by no means determined to be accelerating forever, 
in obvious contrast to standard ACDM cosmology. 

However, in order to relate these results to the observed accelerated expansion of the 
universe, it is crucial to take into account the non-gravitational sector of the theory. 
While Sonner and Townsend [H] considered models comprising only the graviton, the 
axion and the dilaton fields, the purpose of the present paper is to investigate how the 
contribution of cosmological perfect fluids - (dark) matter and radiation - modifies 
the evolution of the AD system. Avoiding any attempt of constructing a fully realistic 
model, we assume that the gravity sector fields couple in a universal metric way to the 
background fiuid. This choice is inspired by the aforementioned ST theories. 

The main result of our analysis is that recurrent acceleration is no longer a generic 
feature of the (modified) AD dynamical system, once we impose the pre-existence of a 
radiation (RDE) and a matter dominated (MDE) era and take into account the finite 
contribution of (dark) matter to the present energy density. 

Our paper is organized as follows. In section [2] we present the cosmological evo- 
lution equations of AD gravity coupled to matter in a universal metric way. Due to 
this coupling any matter field experiences gravitational interactions through the same 
metric Qf^i^, which is conformally related to the Einstein frame metric by a dilaton- 
dependent function. As it will be made explicit in section [2], the dilatonic part of the 
interaction is parametrized by a function Q. The simplest case of minimal coupling 
corresponds to Q = and will be discussed separately in sections [3] and HI the more 
general case of constant Q > is treated in section [5l In section [6l which is devoted to 
our conclusions, we also briefiy comment on the viability of AD models as candidates 
of dynamical dark energy. In the appendices we summarize some basic facts concerning 
dynamical systems terminology and give details of the calculations. 
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2 Axion— dilaton cosmology 



The class of models we are interested in is described by the following action 



•S* — 5'ad + 'S'fiuid 



(1) 



where 




and 



In Eq.([T]) R is the Ricci scalar constructed from the Einstein frame metric g^y, ^{^) 
is an arbitrary function of the dilaton field (to be specified below) and 7, as well as 
A > 0, are real constant parameters. The background fluid sector is described by the 
action S?^aiA■ Hereafter "background fluid sector" refers to the sector of the theory that 
includes all the flelds ^ of the Standard Model (or of one of its possible extensions) 
which we assume to be coupled to the gravitational sector (Sad) by means of the same 
metric g^u = A^($) g^,l,. 

Strictly speaking, by this choice we neglect any non-universal couplings of the AD 
system , in particular interactions with the fleld strength of some gauge fleld, i.e. 
~ F^^F^^ in the case of the dilaton, and ~ a F^^F^^'^ in the case of the axion. In 
fact, such couplings are generically present in any theory which couples AD gravity to 
a matter and gauge sector (see for example fT6]). However, a proper treatment of 
these interactions and their consequences in a cosmological framework is beyond the 
scope of this paper. On the other hand, our analysis applies - more generally - to any 
complex scalar fleld with modular invariant kinetic term. 

The introduction of an exponential potential for the dilaton explicitly breaks the 
S'L(2,M) invariance of the usual AD gravity. Such a potential with 2A = 7 emerges, 
for instance, from a truncation of the Freedman-Schwarz Z) = 4 supergravity theory 
[13] (see also section 4.3 of for an alternative motivation). 

In a flat FRW Universe, 



ds^ = -dt^ + a^{t)dx^ , 
the cosmological equations obtained from the action ([T]) are 



(2) 



^Invariant under 5'L(2, M)— transformations r 



ar-^b 
cr-\-d 



with ad — he — 1 where r = -f ie 2 
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a = -(3//-7<i>)(T, 
1 



H 



(1 + Wfluid)pfluid + + 6-^*^2 



Pfiuid + ^($' + e-^*d^) + e-^* 



where Q{^) is given by 



din A 



(3) 
(4) 
(5) 

(6) 
(7) 



The background fluid energy density pfluid obeys the continuity equation 



Pfiuid 



-3(1 + Wdnid) H-{1- 3Wfl,id) Q($) $ 



Pfiuid, 



with equation of state parameter w^mA ^ [0, |], the limit values corresponding to pure 
matter and pure radiation respectively. 

Introducing the following dynamical variables, 



$2 



X2 



Prad 

3IP 



(9) 



the system (E])^® can be rewritten in autonomous form, straightforwardly generalizing 
the set-up of [T7], 



+ - 1) + [-1x1 + w + Q(i - xl 



2222 
X2 - y - 2; 



dy 
dN 
dz 

In 



■y{xl + xl-y^ + -z^ + l)- 
■z{xl + xl-y^ + l{z^ -1)). 



7x1X2, 



Xxiy, 



(10) 
(11) 
(12) 
(13) 



where N = In a, and Q is assumed to be a positive real constant, corresponding to 
A($) ~ exp(-Q$). 

The given system of evolution equations defines a three-parameter family of dynam- 
ical models with four-dimensional compact phase-space. 



xl + xj + y"^ + z^ < 1, 
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and can, furthermore, be restricted to {x2, y, z) > 0, since the system is invariant under 
change of sign in any of these variables. 

Using the new set of variables, the effective equation of state parameter can be 
conveniently expressed as 

_ Pad + Pfluid 2 , 2 2 , 1 2 

Wcs = ■ = x^+x^-y +-Z . 

Pad + Pfluid <J 

The equation Wes{xi,X2, y,z) = —1/3 defines the boundary of the domain of accelerated 

expansion in phase-space. 

The purpose of the succeeding section is to reproduce the findings of [H] within the 

chosen framework, corresponding to the theory of Eq.([T]) truncated by pfluid = 0. 

3 Axion-dilaton dynamics (/Sfluid = 0) 

Provided a fiat FRW universe and absence of a cosmological perfect fluid, the phase- 
space of the AD system is two-dimensional. We choose it to be spanned by {xi,X2}. 
The Friedmann constraint equation now reads 

2 I 2 I 2 -1 

Xi + x^ + y =1, 
and we can eliminate y from the system: 

'^•^l o__/_2| 2 IN, /^r 2 



— = 3xi{xi + xi-l) + )^-[-lxi + X{l-xi- xi)] , (14) 

dX2 ^ /2 2 



— = 3X2{X^+X2-1) + ^--fXiX2. (15) 

The equation of state is then given by 

Wes = Wad = '2{xl + x^) - 1. 

Each model is characterized by a number of stationary solutions, or critical points, 
of the corresponding autonomous system. (We refer to the appendix concerning a 
brief summary of relevant terminology.) We find the following set of stationary points 
= {xi^s,X2,s), given as functions of the parameters: 



B^,B2: (±1,0), 



G : 




( V6 / A(7 + A)-6 \ 
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fixed point 


existence 


stability 


w 


Bi 


V(7,A) 


stable: 7 < A A > ^6 


1 


B2 


V(7,A) 


saddle point: 7 > 


1 


G 


A < V6 


stable: A(A + 7) < 6 


-i + f 


J 


7 > A A(A + 7) > 6 


stable: 7 > A A(A + 7) > 6 


A+7 



Table 1: Properties of the fixed points of the reduced dynamical system. 



^^3/ i,j e{l,2} 



Properties of the fixed points are displayed in table [H The existence condition can 
be expressed as follows, 

with Xi,X2 real. The stability of a critical point is determined by the eigenvalues of 
the Jacobi matrix 

M 

of the vector function 

F{xi,X2) = l^3xi{xl + xl- 1) + '^^{-jxl + A(l - xl- xl), 

evaluated at Xg. (See the appendix for more details.) We find the following eigenvalues: 

Bi: J ^7, 6-V6A, 



B2: 



G: ^(A(A + 7)-6), 



2(7 + A) I -7±'^72 + 87(7 + A)-^7A(7 + A)2 



The fixed point J is a spiral focus if 

37 (97 + 8A) - 4 7A(7 + A)2 < 0. (17) 
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x[l] 

Figure 1: Model with parameters (A, 7) = (|, l) .Shaded (yellow) area corresponds to 
accelerated expansion. Trajectories are plotted in the plane (xi,X2). 

3.1 Recurrent acceleration 

We now discuss under which conditions recurrent periods of acceleration can be realized 
within the family of dynamical models given by Eqs. (|T^ - fllSp . 

Let us first observe that, according to table [H accelerated expansion is possible at 
the fixed point G - if A < \pl - or at the fixed point J - if 7 > 2A. Then, restricting 
ourselves to the case 7 > 0, we can distinguish three possibilities to realize a model 
which generically allows for periods of accelerated expansion. We give examples of 
phase portraits of the different cases below. 

As we will see, in agreement with |13] , recurrent periods of acceleration are efficiently 
produced by models corresponding to a subset of parameter space where the fixed point 
J is stable and a spiral focus (figure E]). 

3.1.1 G stable 

If A < v^, the attractor G is situated within the domain of accelerated expansion. In 
this case, once acceleration has set in, it will last forever . If \f2 < A < v^, it can be 
a transient phenomenon along a subset of trajectories (figure [T]). 

3.1.2 J stable, G saddle point 

If J is the attractor, the domain of acceleration in parameter space is bounded by 
7 = 2A. The phase-portrait of the system depends crucially on the progress of the spe- 
cial trajectory connecting the saddle point G with the attractor. We will hereafter call 
it the connecting trajectory. If A < V2, both G and J are situated within the domain 
of acceleration in phase-space, and hence the connecting trajectory is completely con- 
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-1.0 -OJ D.D 0.5 1.0 

X[l] 

Figure 2: Model with parameters (A, 7) = (v^, 4). 




-1.0 -D.5 D.D D.5 



Figure 3: Model with parameters (A, 7) = (2,4). 

tained within this domain as well. Any trajectory approaching the connecting one will 
therefore remain inside the acceleration domain once having entered it (see figure [2]). 
If, on the other hand, G is situated outside, recurrent acceleration can be generically 
reahzed, «/the spiral focus J is located close enough to the acceleration boundary, such 
that any trajectory approaching the attractor crosses the boundary repeatedly, as does 
the connecting trajectory (figure [H]). 
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-l.O -0.5 0.0 0.5 L.D 

X[l] 

Figure 4: Model with parameters (A, 7) = (3,4). Only two trajectories are shown. 

L.O — 
0.75 — 

x[2] 0.5 1 

0.25 — 

0.0 — f 1 1 1 1 1 1 1 1 1 1 I I I I I I I 1 1 1 

-L.O -0.5 0.0 0.5 L.D 

x[l] 

Figure 5: Model with parameters (A, 7) = (4,20). Only a single trajectory is shown. 
3.1.3 J spiral focus, G non-existing 

In this case, recurrent acceleration is most generically realized. At A = v^, the fixed 
point G merges with Bi. The dynamical evolution of the system is now totally de- 
termined by the saddle points Bi and -B2, situated at the phase-space boundary, and 
the spiral focus J. (The condition fll7l) is trivially fulfilled in this part of parameter 
space.) Each trajectory winds around the attractor several times, undergoing sub- 
sequent stages of accelerated and decelerated expansion. As figures H] and [5] show, 
this feature is almost independent of the position of the attractor with respect to the 
acceleration boundary. 
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4 Axion-dilaton dynamics in presence of 
cosmological fluids: the case Q=0 

We are now prepared to discuss stationary solutions of the full system ( TT0|) -( |T3i) . but 
setting Q = 0. We find the following set of critical points Xg = {xi^s,X2,s,ys, Zg)'- 



A 

Bi, B2 
C 



E : 



F : 



G ■ 



J : 



(0,0,0,0), 

(±1,0,0,0), 

(0,0,0,1), 

2 

.0, 



^ 0,^,0 



1 - 



A2 



V2A' ' y2A 
A 



-,0 , , 



^/Q /A(7 + A)-6 



7 + A' 



(7 + A)^ 



7 



7 + A' 



The density parameter of the AD system is given by 

f^AD = xl + xl + y. 
Furthermore, we note that now, in the general case, 

^ / 2 , 2 2\ I 2,2 2,^2 

- — (xi + X2-y ) = Wad t u^es = x^ + x^ - y + i:Z . 

"AD <J 

The existence condition reads 

xl + xl + y"^ + < 1, 

with xi,X2,y,z real. The eigenvalues of the Jacobi matrix determining stability of the 
different fixed points are given in the appendix. We display properties of the fixed 
points in table [21 

We find a radiation dominated repeller C, a matter dominated saddle point A, and 
three different AD dominated regimes, B12, G, and J, already present in the reduced 
system (S'fluid = 0). In addition, there are two different scaling solution^. 



See [IHl for a definition. 
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fixed point 


existence 


stability 






A 


V(7, A) 


saddle point 








Bi, B2 


V(7, A) 


unstable 


1 


1 


C 


V(7,A) 


unstable 





1 

Q 


E 


A > 2 


saddle point: 










A > max{2, 27} 


4 

A2 


1 

3 


F 


A > V3 


stable: 










A > max{-\/3, 7} 


3 

A2 





G 


A < V6 


stable: A < 










A A(A + 7) < 6 


1 




J 


A(A + 7) > 6 


stable: 








A 7 > 


7 > A 


1 


A-7 
A+7 



Table 2: Properties of the fixed points of system (JT0l) - (fT3l) . with (5 = 0. 



E, where the energy density of the AD system scales like radiation, and F, where it 
behaves like matter. These two fixed points are characterized by ^ad < 1- 

The subset of fixed points exhibiting X2 = is identical to the set of fixed points 
characterizing single- field models with exponential potential (see [2], and references 
therein). These fixed points correspond to trivial solutions of the axion equation of 
motion (jlj) and our analysis shows that such configurations are indeed stable in a wide 
range of parameter space. 

The existence of stationary solutions with X2 7^ is related to the sign of the fric- 
tion term in the axion equation of motion, i.e. the quantity 3// — 7$. As long as 
"iH — 7$ > 0, the axion evolves toward a configuration where d = and there- 
fore X2 = 0. On the other hand, stability of the fixed point J with X2 7^ implies 
a^ilj > 7~"'"a/3/2, which is equivalent to 3if — 7$ < 0. 

In figures E] and [7] we show the different domains of stability in parameter space of 
both theories, with and without a cosmological fluid background. Most relevant is the 
appearance of the new fixed point F, which is either a stable focus or a saddle point 
in a significant range of parameter space. 

4.1 Recurrent acceleration in presence of background fluids? 

In this subsection we investigate to which extent a perfect fluid background affects the 
possibility of recurrent acceleration. 

Following our discussion presented in section [HI we assume J to be a spiral focus, 
located close enough to the acceleration boundary in phase-space. In other words, we 
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1 2 3 4 5 



Figure 6: Parameter space of the Sonner-Townscnd family of models. Horizontal axis: 
A, vertical axis: 7. In the gray region the fixed point G is stable. The dashed line is 
the existence boundary of G. In the shaded region acceleration is impossible at the 
attr actor. 




1 2 3 4 5 



Figure 7: As above, now the general case. The existence of the new fixed point F is 
indicated for A > -\/3, in the dark gray region F is the attractor. 
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restrict ourselves to a subclass of models, corresponding to the neighborhood of the 
line 7 = 2 A in parameter space. We have to discuss three different cases in turn. 

• A < G saddle point, accelerated expansion at G\ 

• \/2 < A < v^: G saddle point, deceleration at G\ 

• A > v^: F saddle point. 

We will focus on the behavior of the trajectory connecting the saddle point [F or G) 
with the attractor J. Generalizing our previous definition, hereafter we will refer to 
such a trajectory as connecting trajectory. 

Starting with the first case, we note that w^^iCT) > WcsiJ)- This follows from the 
existence condition of J, which can be rewritten as A^/3 > 2A(A + 7)^^. The connecting 
trajectory is entirely contained not only within the domain of accelerated expansion, 
but also in the section of the phase space boundary defined by ^ad = 1- Therefore, in 
this case we find no qualitative difference with respect to section [3l (See figure [2] for 
comparison.) 

In the second case, the limiting trajectory itself crosses the acceleration boundary 
while spiraling onto the attractor. (See figure [3] for comparison.) In particular, any 
trajectory which enters the acceleration domain before approaching the saddle point 
G will experience at least two distinct stages of accelerated expansion. The first stage 
corresponds to the well-known freezing regime [2] of single-field models of dark energy: 
the dilaton field velocity remains close to zero due to the Hubble friction term dominat- 
ing the equation of motion. The second stage is reached, when the trajectory re-enters 
the domain of acceleration in approaching the connecting trajectory and the regime of 
the late time attractor. If J is located inside the domain, accelerated expansion will 
continue forever. 

In the third case the saddle point F is dynamically relevant. Moreover, at F 
Wcs = Wad = and the expansion is either dominated by matter or by the AD sys- 
tem scaling like matter, with Qad < 1- Due to the finite contribution of matter to the 
total energy density of the universe, we have now - in contrast to the previous cases - 
Wes > Wad along the connecting trajectory. Thus, even if the connecting trajectory 
oscillates around wad{J) ~ — 1/3 before reaching the spiral focus, recurrent accelera- 
tion is not implied. In particular, in any model with Wes{J) > it will never enter 
the domain of acceleration at all. However, we cannot exclude the possibility of Wes 
crossing the acceleration boundary more than once, if at the attractor Wes{J) < —1/3. 



13 



4.2 Numerical examples 

The conclusions we drew in the previous subsection can be circumvented by resorting 
to a very special choice of initial conditions. For instance, the scenario of [Hj can be 
recovered by setting 

^in = 0, {xj + xl + = 1. 

However, since we are ultimately interested in models which are able to reproduce 
qualitatively the standard evolution of the universe, as it can be reconstructed from 
cosmological observations 0, we will hereafter only consider trajectories which are (at 
least marginally) consistent with the concordance ACDM cosmology. Furthermore, we 
will henceforth identify I^ad with Qde and wad with w^e, where the subscript DE 
refers to dark energy. 

It is a well-known, serious problem of dynamical dark energy models that their late 
time evolution typically still depends on initial conditions: In single-field models, for 
instance, the energy scale of the potential, i.e. in our notation ?/in, has to be fine tuned 
to satisfy ^DE.today ~ 0.75. In the present case one has to deal with an additional 
sensitivity on X2^in, which will be explained below. 

We will discuss these issues on the basis of two numerical examples, which are both 
characterized by an attractor solution given by fixed point J, preceded by a saddle 
point, which is G in the first case, and in the second one F. We have specified the 
respective trajectory by imposing initial conditions for the AD system at a temper- 
ature of 0{1) MeV when (f2i.ad/^mat)in ~ 10^, thereby ensuring the validity of our 
classical description. We assume the scalar fields to have already reached the freezing 
regime, relying on the fact that scalar field kinetic energy scales as a~^. This leaves the 
possibihty of a stage of kination (TB] during a preceding epoch of higher temperature. 

As it can be seen in figures E] and [H we find, in both cases, three successive evolu- 
tionary stages of the AD system: The first one, the freezing regime, is associated with 
the radiation dominated epoch (RDE) if F is the saddle point, or continues during the 
matter dominated epoch (MDE), if G is the saddle point. Thereafter the system enters 
the regime of the saddle point, which in both cases lasts for a significant number of 
e-foldings of expansion. Finally, there is the late time attractor regime. Remarkably, 
our present situation corresponds to the transition between stage one and two in one 
case (saddle point G), and two and three in the other (saddle point F). This is re- 
lated to the significant discrepancy in ?/in in our two examples (see the figure captions). 
Changes in i/i^ affect in particular the termination of the freezing regime, while the 
choice of X2,m determines the duration of the saddle point regime. 

The value of X2 decreases monotonically along a given trajectory as long as Xi < 
1^3/27""^ and turns to increase when xi > \Jdt/2'^~^ (which is already true at the 
saddle point F). The saddle point regime ceases, once X2 has increased sufficiently to 
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perturb Xi away from its fixed point value: A non-zero axion kinetic term contributes 
to the effective potential in the dilaton equation of motion, Eq. (I3l). 

On the contrary, the initial condition xi^in influences just the early stage of the 
dynamical evolution: A set of trajectories differing only in Xi^in first converge toward 
the saddle point before they start to approach the attractor, thereby washing out any 
dependence of the late time evolution on xi^n- The corresponding evolutionary path in 
phase-space is therefore completely determined by the connecting trajectory, provided 
we can safely assume X2^in to be sufficiently small. 

To demonstrate the impact of a cosmological background fluid on the AD dynamics, 
in particular concerning the phenomenon of recurrent acceleration present in the case 
5'fluid = 0, we have also plotted - for comparison - the evolution of the system's 
equation of state in absence of the fluid background (observe the green lines in the 
figures). The admixture of a perfect fluid component to the initial composition has 
two effects on the evolution of the scalar fields: First, the reduction of f^DE.m by a huge 
factor (9(10~^°) allows for the existence of a freezing regime during RDE and/or MDE. 
Secondly (see figure [H]), due to the existence of the scaling saddle point F, the first few, 
large amphtude oscillations of the equation of state are suppressed and partly replaced 
by oscillations around the saddle point value, leaving only rapid, small amplitude 
oscillations around the attractor value. Needless to say, this kind of oscillations in 
the DE and effective equation of state are - at low redshift - already disfavored by 
observational data [2], [19] . 

We have to conclude that the characteristic feature of recurrent acceleration, as illus- 
trated in figure [5], disappears if we allow for a perfect fluid contribution dominating the 
earlier stages of evolution. In particular, the existence of a scaling solution preceding 
the spiral focus regime reduces number, amphtude and period of possible oscillations 
in Wes crossing the acceleration boundary. 
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Figure 8: Model with parameters (A, 7) = (|,4). Evolution of wde (blue, dark curve), 
Weff (yellow, light curve), with N := — ln(l + z), N = referring to the present. The 
trajectory is specified by initial conditions (xi, X2, y)in = (10~^^, 10~^^, 10"^''). For com- 
parison, the green curve shows the evolution along a trajectory with fimat = ^rad = 0, 
as in the Sonner Townsend case, but the same ratio between initial kinetic and potential 
energy of the scalar fields. 



16 



1.0 
0.8 
0.6 
0.4 
0.2 

W 0.0 
-0.2 
-0.4 
-0.6 
-0.8 
-1.0 

-20 -15 -10 -5 5 

N 

Figure 9: Model with parameters (A, 7) = (4, 20). Evolution of i^de (blue, dark curve), 
Weff (yellow, light curve), with := — ln(l + z), N = referring to the present. The 
trajectory is specified by initial conditions (xi, X2, y)m = (5 x 10~^^, 5 x 10"^*^, 5 x 10~^°). 
For comparison, the green curve shows the evolution along a trajectory with flumt = 
^rad = 0, as in the Sonner Townsend case, but the same ratio between initial kinetic 
and potential energy in the scalar field sector. 
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5 Axion-dilaton dynamics in presence of 
cosmological fluids: the case Q > 

The set of equations ( fT0i) -( |T3l) . with Q allowed to be non-zero, defines a three-parameter 
family of models, each characterized by a set of fixed points in a four- dimensional 
compact phase-space. These fixed points are: 

A: |^y|g, 0,0,0 

Bi,B2: (±1,0,0,0) 
C: (0,0,0,1), 

D: ^ ' 
E : 




^.o„/H^?^g:ll±l.o 



A-Q' 'V 2(A-Q)2 



G: l^-0,;i-^,0 



7 + g'V 2(7 + g)2 



J 



/A(7 + A)-6 



7 + A'V (7 + A)^ 




With respect to the case Q = 0, we find two additional stationary points, namely 
D and H. The fixed point D is associated to RDE. Properties of the fixed points are 
displayed in table El 

In figures [TO] and [TT] we show two sections of parameter space, with Q = 1/2 and 
Q = 1 respectively, to cover the different possibilities of stable fixed points. Due to 
the positivity of Q, the AD energy density gets enhanced at expense of the matter 
sector. If Q is sufficiently large, not a single fixed point remains with Wes equal or at 
least close to zero, indicating suppression of MDE. However, we have to note that such 
large values of Q are unphysical because of the existing bounds on a universal metric 
coupling between matter and gravity [SI [20] • 

Moreover, let us emphasize that increasing Q does not re-establish recurrent acceler- 
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f.p. 


existence 


stability 






A 




stable: 
















B, 




saddle point: 








V(7, A, Q) 


A>V6AQ>y^A7<0 


1 


1 


B2 




unstable 






C 


V(7,A,g) 


unstable 





1 
3 


n 






1 


1 
3 


E 


A > 2 


stable: 27 < A < A.Q 


4 
A2 


1 
3 


F 




stable: Q < ^ 










A A > max{4(5, 2(5 + 7} 


(Q-A)2 


Q 
\-Q 


G 


A < V6 


stable: 










A^ < mm{4, 3 + QA, 6 - 7A} 


1 


-i + ¥ 


H 


7 > max{0, ^ - <5} 


stable: A > 7 + 2(5 








A 7 > 2Q 


Q 

7+Q 


Q 

7+Q 


J 


A(A + 7) > 6 


stable: 








A 7 > 


A < mm{27, 7 + 2(5} 


1 


A-7 
A+7 



Table 3: Properties of the fixed points in the case (5 > 0. 



ation. As in the case (5 = 0, generic trajectories which converge toward the spiral focus 
J approaching the connecting trajectory, always spend a certain number of e- foldings 
close to a saddle point, where either radiation or matter dominate. Depending on the 
parameter values, the relevant saddle point is either F, A, H, D or E. 
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1 2 3 4 5 



Figure 10: Parameter-space in the case Q = ^. Horizontal axis: A, vertical axis: 
7. Regions of stability of the various fixed points are indicated. The full line is the 
stability boundary of J. F can only be stable if2<A<7/2. For A > 7/2 we find new 
attractors A (if 7 < 5/2) and H. In the shaded region acceleration is impossible at the 
attr actor. 




12 3 4 



Figure 11: As above, now the case Q — 1. Attractors F and A are replaced by E (if 
2 < A < 4) and D, both corresponding to radiation era. 
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6 Discussion and conclusions 



We have investigated a family of flat FRW cosmological models in D = A, focusing on 
the dynamics of a scalar (dilaton) and a pseudo-scalar (axion) partner of the metric field 
in presence of cosmological background fluids. Neglecting any specific interaction term 
of the axion and the dilaton, which are present in particle physics or string motivated 
models, we have considered the - simplified - scenario of universal metric coupling. 
Motivated by recent results of [H] , where recurrent acceleration was verified in a large 
class of AD models in absence of matter or radiation, we have analyzed the possibility 
to combine the phenomenon of recurrent acceleration with a cosmological background 
evolution in concordance with ACDM cosmology. In particular, we have considered 
the viability of the AD system as a candidate of dynamical dark energy. 

Using a different choice of dynamical variables, well-suited to the more general case 
with background fluids, we were able to reproduce the findings of [H] (corresponding 
to Sfiuid = 0) . The new feature of models with Sfiuid 7^ is the existence of a stationary 
cosmological scaling solution within a wide range of parameter space. As long as Q is 
(close to) zero, the relevant fixed point is F, corresponding to matter dominance. In 
particular, if being a saddle point, F influences the evolution along generic trajectories 
in such a way that the phenomenon of recurrent acceleration, as observed in jTlj, is 
reduced to small amplitude oscillations of the equation of state at low redshift. 

We can certainly realize a model (by specifying parameters), which is able to re- 
produce the present stage of cosmic evolution as being a transient phenomenon. This 
is already possible in the single-field case, corresponding to the constant axion sce- 
nario within our dynamical system: The dilaton remains frozen, due to Hubble friction 
dominance, right up to the present, and later-on evolves toward an attractor solution 
exhibiting > —1/3. On the other hand, if the axion dynamics is non-trivial, it 
is also possible to get a second accelerating stage in the future, which will then be 
ever-lasting. 

Though we have to conclude that recurrent acceleration as described in [T3] is not 
relevant to dark energy model building, we have discovered a different, interesting 
possibility instead, which we may call the F J scenario. Does it provide a viable 
model of dynamical dark energy? 

The co-existence of a matter-dominated scaling solution and a dark energy domi- 
nated accelerating solution is generally considered to be a very attractive feature of 
dynamical DE models [21 [IE]- In the single-field case [17], the existence of F is excluded 
by the stability condition of the fixed point G, which is the only available candidate 
to achieve late time acceleration. However, in presence of an axion field, the scaling 
regime of F can be succeeded by a stage of accelerated expansion, represented by the 
fixed point J. During radiation and matter dominated epochs, the evolution of the 
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AD sector is determined by the saddle point, while our present situation corresponds 
to the transition toward the late time attractor solution. Unfortunately, the onset of 
this transition is subject to a certain fine tuning of initial conditions. 

One single fine tuning of the scalar potential energy scale is always mandatory in 
dynamical models of DE, since it corresponds to setting the cosmic clock. The re- 
quirement of naturalness [21] strongly disfavors models which need a fine tuning of 
the same order as in the cosmological constant case. In this respect, two-field models 
incorporating the F J scenario are certainly promising: In our numerical example, 
the potential scale is initially set to 0(10^^) x J^aIimcV- Even larger values are possible, 
but have to be compensated by reducing X2,in, in order to keep the cosmic clock tuned. 
We need to emphasize, however, that extremely small values of the axion field velocity 
are by no means unphysical. In fact, the most natural solution is a constant axion, as 
long as the Hubble rate is large enough to keep the friction term positive. 

To conclude, let us stress an intriguing feature of the class of models under consid- 
eration. Single-field potentials are typically required to be extremely flat in order to 
satisfy observational constraints. However, in the F ^ J scenario this is not true: 
there is no upper bound on the parameter A determining the potential slope. This 
aspect opens up new possibilities for dark energy model building which we believe to 
be promising enough to motivate further investigation. 
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A Dynamical systems terminology 



We consider a system of n first order ordinary differential equations (ODE), 

Xi^ fi{xu...,Xn), (18) 

wliicfi is called autonomous if none of the n functions fi explicitly depends on time. A 
solution of the system is given in terms of a trajectory in phase space, 

t ^ X(t) := (xi(t),...,Xn(t)), 

determined by choice of initial conditions X{tinit)- 

A point Xs :— {xi^s, Xn,s) is said to be a critical, stationary or fixed point if 

MX,)=0 yi<n, 

and an attractor if there exists a neighborhood of the fixed point such that every 
trajectory entering this neighborhood satisfies the following condition: 

lim X{t) = Xs. 
Now consider small perturbations around the critical point, 

Xi — X-i^s ~t~ ^Xi- 

Linearizing the evolution equations we obtain a system of first order ODE linear in the 
perturbations, 

—Sxi = ^Mij5xj, (19) 
j 

where 

The general solution of this system is given by 



x=x. 



= J^Qfee'^'^*, 

fc=i 

where Cik are integration constants and /Xfc the eigenvalues of the Jacobi or stability 
matrix M, which we have assumed to be distinct for simplicity. Obviously the pertur- 
bation will decay if each Hk has negative real part. 

The critical points of a dynamical system can be classified in terms of the eigenvalues 
of the corresponding stability matrix. An attractor is characterized by the requirement 

Re[fik] < y k < n, 

and called spiral focus if at least one pair of eigenvalues is complex and stable node else. 
Furthermore we will use the terminus saddle point if and only if M has one eigenvalue 
with positive real part. In any other case we call the fixed point unstable. 
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B Jacobi matrix eigenvalues 



In the following we list the Jacobi matrix eigenvalues at the different fixed points 
the models with background fluids. 



B.l Case Q = 



3 3 3 1 
A : — , — . -. — , 

2' 2' 2' 2' 



Bi: 3, 1, W^7, 



B,: 3, 1, -J^7, 3+ J^A, 



C: 2,-1,-1,1, 



^ _1 3(7 - A) 3 ^ V24A2 - 7X^ 



2' 2A ' 4 V A2 
G : ^(A^ - 6), i(A(A + 7) - 6), A^ - 3, ^(A^ - 4), 



J : 3 



37 
7 + A' 



2(7 + A) I ' V ' ' 3 
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B.2 Case Q ^ 



B,: 1, 3-V6Q, y|7, 3 - ^^A, 



B^ : 1, 3 + V6Q, -W^7, 3 + 



C: 2,-1,-1,1, 



S : 1 



2g' 2 g2 

4Q _ 27 1 / V64A^ - 15An 
A' ^A'2lv A2 

A-4g 3 / ^ , 7 + Q 



2(A-g)' 2 V A-gy ' 



-(A - 2g)(A - g) ± W (A - g)2[24 - ta^ - i2Ag + 2og2 + ^Ag(A - qy 



4(A-g)2 \ ^ ' V ' ^ ^ ' ^ ' 2, 

G : \{\'- 6), i(A(A + 7) - 6), - 3 + A(A - g), \{\' - 4), 



TT . 1 37 3 A-g 

2(7 + g)' 2^' 7 + g^' 

^ (-37 ± - 247g[2(7 + Qf - 3]) , 



4(7 + Q) 



J: 3(l-^il±^Vl- 



7 + A / 7 + A 



-7 ± \/72 + 87(7 + A) - ^7A(7 + A)2 



2(7 + A) V V ' ^ 3 
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